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Abstract 

We consider in parallel three one-dimensional spin models with kinetic con- 
straints: the paramagnetic constrained Ising chain, the ferromagnetic Ising chain 
with constrained Glauber dynamics, and the same chain with constrained Kawasaki 
dynamics. At zero temperature the dynamics of these models is fully irreversible, 
leading to an exponentially large number of blocked states. Using a mapping of 
these spin systems onto sequential adsorption models of, respectively, monomers, 
dimers, and hollow trimers, we present exact results on the statistics of blocked 
states. We determine the distribution of their energy or magnetization, and in 
particular the large-deviation function describing its exponentially small tails. The 
spin and energy correlation functions are also determined. The comparison with 
an approach based on a priori statistics reveals systematic discrepancies with the 
Edwards hypothesis, concerning in particular the fall-off of correlations. 
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1 Introduction 



In a great variety of systems, such as structural glasses, spin glasses, and granular ma- 
terials, the dynamics at low temperature or high density is so slow that the system falls 
out of equilibrium [|I]. For long, glassy dynamics has been described as a slow motion 
in a complex energy (or free energy) landscape f^, with many valleys separated by bar- 
riers. Several approaches have been proposed in order to make this heuristic picture 
more precise. Valleys thus appear under various names (and with various definitions) 
in different contexts: metastable states, TAP states 0, pure states inherent struc- 
tures 0, quasi-states 0. From a dynamical viewpoint these concepts are not equiva- 
lent [0]. Metastability is indeed unambiguously defined for mean-field models only, where 
metastable states have infinite lifetimes, as barrier heights diverge with the system size. 
For finite-dimensional systems with short-range interactions, barrier heights and valley 
lifetimes are always finite at finite temperature, so that metastability becomes a matter 
of time scales p. 

Once these valleys are appropriately defined, one can estimate their number at fixed 
energy (or free energy) density E. This number generically grows exponentially with the 
system size, as 

Ar(iV) ~exp(iV5ap(S)), (1.1) 

where S^piE) is the configurational entropy or complexity. The subscript 'ap' (a priori) 
refers to the fact that, when counting valleys, each of them appears with the same weight. 
In this a priori ensemble, valleys are combinatorially equivalent. 

In contrast, a key question concerns the dynamical weight of each individual valley. 
Do all the valleys play a similar role in the dynamics? This question arises for instance 
when a system is instantaneously quenched into the glassy phase, starting from a disor- 
dered configuration. Assuming this initial configuration is chosen at random, does the 
system sample all the possible valleys with a given final energy E with equal statistical 
weights, i.e., with a uniform or fiat measure, or, to the contrary, does the size of the 
attraction basin of each valley really matter? The same question is also relevant in an- 
other situation commonly referred to as tapping ||, |1^, |1^. Under tapping a granular 



material continuously jumps from a blocked configuration to a nearby one. Does the 
non-equilibrium steady state thus obtained admit a statistical description in terms of a 
fiat ensemble of blocked configurations? 

The answer to this question for the first situation (relaxational or aging dynamics) 
is positive at least in some mean-field models, where valleys are known to be explored 
with a fiat measure [^. The concept of ergodicity, and the resulting thermodynami- 
cal construction, therefore hold, as in equilibrium situations, up to the replacement of 
configurations by valleys. The configurational temperature Tap, defined by 

^ = (1.2) 

Tap dE ' ^ ^ 
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has a thermodynamical meaning. It also coincides with the effective temperature involved 
in the generalized ffuctuation-dissipation formula in the appropriate temporal regime. 

Besides the mean-field geometry, another physical situation where metastable states 
are unambiguously defined is the zero-temperature limit, where no barrier can be crossed 
at all. Valleys are just blocked configurations under the chosen dynamics. For instance, 
for an Ising model with single-spin dynamics, a valley is a configuration where each spin 
is aligned with its local field, provided the latter is non-zero. 

In the context of granular materials, Edwards |]12| proposed to describe the slow com- 
paction dynamics by means of a flat ensemble average over all the blocked configurations 
of the grains with prescribed density. Extending the range of application of this idea, 
the so-called Edwards hypothesis consists in assuming that all the valleys with a given 
energy density are equivalent. This hypothesis has two consequences. First, the value of 
an observable can be obtained by a fiat average over the a priori ensemble, or Edwards 
ensemble, of all those valleys. Second, the temperature Tap of (|1.2D, also known as the 
Edwards temperature, has the usual thermodynamical meaning of a temperature. 

The present paper is devoted to the analytical study of the zero-temperature dynam- 
ics of simple one- dimensional systems with kinetic constraints. Kinetically constrained 
models have been the subject of numerous investigations [|13|, |l^ |15|, |16|, |l^ |18|, |19|, ^ ^ 



22, 23, pi, UB, E™. Here we specialize to Ising chains without frustration nor quenched 



disorder, namely paramagnetic chains (CIC) |T^, |T5|, [T6|, |T^ (section 2) and ferromagnetic 



chains (ffSl EO, 21, |22| (sections 3 and 4). The common feature of all these models is 



the irreversible nature of their zero-temperature dynamics: each spin flips at most once 
in the whole history of the system. 

The central goal of this paper is to obtain exact results on the statistics of the blocked 
configurations reached by these systems, pursuing the efforts made by the authors of some 
recent works [|1^, |19|, ^ |23|, |2^, ^ . Thus doing we are able to critically revisit 

the questions raised above. In the derivation of these results we will take advantage of 
the fact that the zero-temperature dynamics of these models can be rephrased in terms 
of random sequential adsorption (RSA) or cooperative sequential adsorption (CSA) [ p7| , 
for which analytical techniques are available in one dimension. A blocked configuration 
thus appears as a jammed state of the corresponding RSA or CSA model. The Edwards 
hypothesis is, in this specific context, akin to the question raised long ago PB| whether 
RSA configurations are equilibrated or not. Our results confirm that the answer is 
negative. A more extensive discussion will be given in section 5. 



2 Constrained Ising chain 
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2.1 Definition of the model 



Constrained Ising chains (CIC) [|I^, |T^, |r7| are among the simplest examples of 



kinetically constrained models |13|- Although they have trivial equilibrium properties, 
the presence of kinetic constraints leads to slow dynamics at low temperature, and to 
metastability at zero temperature. 

Consider a paramagnetic Ising chain, made of independent spins cr„ = ±1, submitted 
to a positive unit magnetic field. This model has a Hamiltonian 

n = -Y.^n^ (2-1) 

n 

with a unique ordered ground state, where all the spins are up (cr„ = +1). 

Kinetic constraints are introduced as follows. Consider single-spin-flip dynamics with 
rates 

W{an -(Tn) = min(l,e-2'^"") Wo((^n-i, 

The first factor in the right side is the Metropolis acceptance rate, ensuring detailed 
balance with respect to the Hamiltonian 7i at temperature T = 1/ (3. The second factor 
imposes a kinetic constraint: the flipping rate of a spin o"„ depends on its environment, 
i.e., on the value of its neighbors o"„_i and (Xn+i- Let us make the choice |23 



Wo{an-i, (Jn+i) = avn-i + (1 " a)r„+i, (2.2) 

with the notation 



I- an 

Tr, = 



2 

The parameter < a < 1 allows to interpolate between known limiting cases. For 
a = 1/2, the constraint factor is (r„_i + Tn+i)/2, i.e., half the number of neighboring 
down spins. The symmetrically constrained chain (SCIC) ^ is thus obtained. For 



a = 0, the constraint factor is r„+i: the spin an can only flip if its right neighbor is down. 
The right asymmetrically constrained chain (ACIC) |T^, is thus recovered. Similarly, 
the left ACIC is obtained for a = 1. 

At zero temperature, the dynamics of the CIC simplifies drastically. An initially up 
spin {an = +1) remains up forever, while a down spin (cr„ = —1) can flip at most once, 
according to stochastic rules depending on its two neighbors: 

> — I — (rate 1), 

h^- + + (rate a), (2.3) 

H > + -\ — (rate 1 — a). 

For a flnite chain of N spins, the dynamics stops after a flnite jamming time Tjv, which 
depends both on the initial conflguration and on the history of the chain. The jamming 
time will be shown to grow as T/v ~ In N, up to finite fiuctuations given by extreme- value 
statistics. The system is thus left after a finite time T/v in a jammed or blocked state. 
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This state, which is an attractor for the dynamics, is a spin configuration where each 
down spin is isolated, i.e., surrounded by two up spins. The blocked configuration thus 
obtained depends on the parameter a, on the initial configuration of the chain, and on 
its whole stochastic history. 

The problem may be equivalently described as an irreversible process of particle ad- 
sorption. Consider indeed down spins as representing empty sites (o), and up spins as 
representing occupied sites (•): 

{a = +1 r = •. 

The zero-temperature dynamics of the CIC thus maps onto a problem of particle adsorp- 
tion, where individual particles (monomers) are irreversibly deposited according to: 

o o ^ o • o (rate 1), 
o • ^ o • • (rate a), 
o o ^ • • o (rate 1 — a). 

The deposition rate at site n depends on the occupation state of both neighboring sites. 
We are thus facing a cooperative sequential adsorption (CSA) model The limit 



(jamming) coverage of this model is related to the mean magnetization per spin M(oo) 
in the blocked configurations by 

P.(.) = i-W» = -^'"<°°' 



2 

In the following, we will use the language of spins, magnetization, spin correlations, and 
so on, leaving the picture of particle deposition for illustrative purposes only. 



2.2 A priori statistics 

We have shown that the attr actors of the zero-temperature dynamics of the CIC, for any 
value of the parameter a, are the spin configurations where each down spin is isolated, 
i.e., surrounded by two up spins. 

A natural statistical description of these attractors is provided by the a priori ensem- 
ble, or Edwards ensemble, as explained in the Introduction, where all the blocked spin 
configurations are taken with equal weights. 

For a finite chain of spins, consider the restricted ensemble of blocked configurations 
for which exactly n spins are down. Their magnetization M is such that NM = N — 2n, 
with < 72 < N/2, hence < M < 1. The number of such configurations reads 

MiV.n)=('^7+'). (2.4) 

Indeed this is the number of ways of inserting n down spins in the N — n + 1 spaces 
made available by the presence of A^ — n up spins, with at most one down spin per 
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space. This number grows exponentially, according to (|1.1| ), where the a priori entropy 
reads [g3], ^ 

^ap(M) = -Mln(2M) + ^ ln(l + M) - ^~^ ln(l-M). (2.5) 

One can also consider the full (or unrestricted) ensemble of all the blocked configura- 
tions, irrespective of their magnetization. The number Af{N) of such configurations can 
be determined as follows. For a chain of > 3 spins, a configuration either ends with 
(H — ) (there are Af{N — 2) such configurations) or with (+) (there are A/'(A^ — 1) such 
configurations). We thus obtain the recursion relation //{N) = //{N — 1) +N'{N — 2), 
with Af{l) = 2, Af{2) = 3, hence 



Af+2, 



where F/v are the Fibonacci numbers. This expression is also equal to the sum of (|2T 
for n ranging from to N/2. It grows as A/'(A^) ~ exp(A^S'*p), with 

= In $ = 0.481212, (2.6) 

where $ = (1 + ^/E)/2 is the golden mean. The result ( p.6|) is the maximum value of the 
function S'ap(M) ( p.5|) . This maximum is reached for 

1 

= = 0.447214, (2.7) 

which is therefore the typical a priori magnetization of a blocked configuration. 

The distribution of the number n of down spins in the a priori ensemble is given by 



Af{N, 



n] 



" M{N) ■ 

For a large sample {N ^ 1), the probability density of the magnetization M is therefore 
given by an exponential estimate of the form 

/(M)~exp(-iVS,p(M)), (2.8) 

with 

S,p(M) = 5:p-5,p(M). (2.9) 

The result ( p.8|) has the form of large-deviation estimates in probability theory, which 
hold e.g. for the arithmetic mean of independent random variables. The large-deviation 
function (or entropy function) Eap(M) will be plotted in Figure ^j. It vanishes quadrati- 
cally near M = M* as 

S,p(M)^c(M-M^)% c=^. 

o 
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The bulk of the a priori distribution of M is therefore asymptotically a narrow Gaussian 
around M*, with a scaled variance given by A^VarM ^ l/(2c) = 4v^/25 = 0.357771. 

The a priori entropy can alternatively be evaluated by the transfer-matrix method [p9| . 
For a finite chain of spins, we introduce the characteristic function of the magnetization 

c 

where the sum runs over all the blocked configurations C. Note that NM{C) is the 
opposite of the total energy of the configuration, according to the Hamiltonian (|2.1|) , 
so that Z]sr{(3) coincides with the usual partition function of the model, at a fictitious 
inverse temperature /?. 

The partition functions of a finite chain of spins, labeled by the prescribed 
value aN = ±1 of the last spin, obey the recursion 

/ 7+ \ / 7+ \ 



where the 2x2 transfer matrix 



r 



' e^ e^' 
e-^ 



has eigenvalues 



2 

The entropy S'ap(M) is then given by a Legendre transform. We have indeed 



Evaluating the integral by the steepest-descent method yields the 'thermodynamical' 
relationships 

lnA+(/3)-Sap(M)=/?M, M = ^^, ^ = (^.lO) 
which yield in the present case 

M = , e^ = 



_ 4 + e^^ VI - M 

and allow to recover ( p. 51) . 

The spin correlation function Cn = (o"o(t„) can also be evaluated in the a priori 
ensemble at fixed magnetization by the transfer-matrix method ||29|. We have, for n > 
in the bulk of an infinitely long chain, 

Cn = = {{L,\S\R,)f + {L,\S\R^){L_\S\R,) (^Y. 
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In this expression, S = diag(+l, —1) is the spin operator, while 

are the left and right eigenvectors of T associated with the eigenvalues X±. We have 
consistently M = (L+|5|i?+). After some algebra we obtain the following expression. 



involving only the magnetization M [ITH 



= C„-M2 = (l-M2)(-^^j . (2.11) 

The connected correlation function thus exhibits an exponential decay, modulated by an 
oscillating sign. 

The full ensemble of blocked configurations is obtained by setting /3 = in the above 
results, which indeed corresponds to taking a fiat average over all blocked configurations. 
This prescription amounts to replacing the magnetization M by its typical value M* ( |2.7| ). 
We thus obtain in particular 

1 4 / 1 \" 
Cr° = C„ - - = - (-^^) . (2.12) 

We end up by mentioning that the blocked spin configurations considered so far are 
the degenerate ground states of the antiferromagnetic Ising chain in a constant magnetic 
field h = 2J > [BDI, whose Hamiltonian reads 



H = (T„cr„,+i - 2 J ^ an- 

n n 

As a consequence, the above expressions are exact results for the latter model at equi- 
librium at zero temperature. This is one of the simplest models with a non-zero entropy 
at zero temperature, given by (|2.6|). 



2.3 Dynamics of cluster densities and magnetization 

We now turn to the exact analysis of the zero-temperature dynamics of the CIC, starting 
with the mean cluster densities and magnetization. 

We consider an uncorrelated magnetized initial state, given by 

f a„(0) = -1, r„(0) = 1 (o) with prob. p, 

\ a„(0) = +1, r„(0) = (•) with prob. 1 - p, ^ ' ^ 

so that the mean initial magnetization reads M(0) = 1 — 2p. 

For p < 1/2, the initial state (|2.13| ) is the equilibrium state of the Ising chain with 
Hamiltonian ( p.l|) at inverse temperature 

1 1 -V 

/3o = -ln ^. (2.14) 

2 p 
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In particular, a random (unmagnetized) initial configuration, i.e., p = 1/2, corresponds 
to infinite temperature, i.e., Pq = 0. 

It is a common feature of one- dimensional RSA and similar problems that the 



densities of certain patterns, including active clusters, obey closed rate equations. Con- 
sider clusters of exactly i > 1 consecutive down spins. Their density per unit length at 
time t reads 

peit) = ((1 - ro)ri . . . r,(l - r.+i))^ = o_^.), (2.15) 

e 

and the mean magnetization of the chain is given by 

M{t) = l-2j2ipe{t). (2.16) 

Because zero-temperature dynamics is fully irreversible, the densities pe{t) obey rate 
equations, which can be derived as follows. Clusters of length i = 1 are inactive. Consider 
a cluster of length i > 2, renumbering its sites as n = 1, . . . , The spin ct„ can fiip from 
down to given by (|2.3| ), thus generating one or two smaller clusters of the 

following length 

' n = 1 (rate 1 — a), one cluster: ii = i — 1, 

2<n<i-l (rate 1), two clusters: £^ = n - 1, £2 = i - n, (2.17) 

n = i (rate a), one cluster: ii = i — 1. 

Gathering the contributions of all these events, we obtain the rate equations 

= _(£ _ l)p,(t) + pe+,{t) + 2 Y: Pk{t) (2.18) 

^'^ k>i+2 



for i > 1, irrespective of the value of the asymmetry parameter a. The initial state ( p.l3|) 
yields p^(0) = (1 — pYp^. 



A simple way of solving the rate equations (|2.18| ) consists in making the Ansatz 

peit) = a{t) z{tY (2.19) 
for i > 1. We obtain successively dz{t)/dt = —z{t), with z{0) = p, hence 

z{t) = pe-\ (2.20) 
and da{t)/dt = a(t)(l + z{t)^)/{l - z{t)), with a(0) = (1 -p)^ hence 

a{t) = e\l - pe'^Y exp{p{e-^ - 1)), 

so that finally 

p^(t) = (1 - pe-^y exp(p(e-* - 1)) pV^^-^)*. (2.21) 
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As expected, only inactive clusters of length i = 1 survive in the final states, and their 
density reads 

Pi(oo) = pe~^. 



Equation (|2.16|) yields 
and especially 
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M{t) = 1 - 2j9exp(p(e-* - 1)), 
M{oo) = 1 - 2pe-P. 



(2.22) 
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Figure 1: Plot of the final energy E{oo) against the initial energy E{0). Full line: CIC (|2.23| ). 
Dashed line: ferromagnetic chain with constrained Glauber dynamics ( ^.71) . Long-dashed line: 



ferromagnetic chain with constrained Kawasaki dynamics (4.5). 



The mean final magnetization of the blocked states reached by the dynamics thus 
depends on the parameter p characterizing the initial state. This non-trivial dependence 
demonstrates that the dynamics is not ergodic. For an initial state close to the ground- 
state {p 0, i.e., M(0) 1), the behavior M(oo) ^ M(0) + 2p'^ is easily explained 
in terms of clusters of two down spins: the density of these clusters scales as p^, and 
only one of the two spins will fiip. For a random (unmagnetized) initial configuration 
{p = 1/2, i.e., M(0) = 0), we have 



M(oo)p=i/2 = 1 - e^^/2 = 0.393469. 



(2.24) 
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As this number is the final magnetization of a typical initial state, it is natural to compare 
it to the prediction ( p.7| ) of the a priori ensemble. This comparison will be presented in 
Table 1. Finally, for the ordered initial state where all spins are down {p = 1, i.e., 
M(0) = —1), we have the smallest possible value of the final magnetization: 

M(oo)p=i = 1 - 2e-^ = 0.264241. (2.25) 

Figure |l| shows a plot of the final energy E{oo) = — M(oo) against the initial one, E{0) = 
— M(0), for the present model, as well as for the ferromagnetic chain with constrained 
Glauber and Kawasaki dynamics (see sections 3 and 4). 

To close up, we present an analysis of the distribution of the jamming time of a 
large but finite system of spins, a question which does not seem to have been considered 
in previous works on RSA. Equation ( p.21| ) shows that the late stages of the dynamics are 



governed by an exponentially small density of surviving clusters made of two down spins, 
P2{t) ~ ae~*, with a = p^e~^. The dynamics can therefore be effectively described by a 
collection of aN such clusters, each cluster decaying exponentially with unit rate, when a 
down spin fiips. The jamming time T^r is the largest of the decay times of those clusters. 



For a large sample, it is therefore distributed according to extreme-value statistics [pT 
Setting 

Tn = HaN) + Xn, (2.26) 

we find that the fiuctuation remains of order unity, and that it is asymptotically 
distributed according to the Gumbel law 

/(X)=exp(-X-e-^). (2.27) 

We have checked this prediction by a numerical simulation. Figure ^ shows a histogram 
of the observed jamming time T/v for 10® samples of A^ = 1000 spins, starting with a 
random initial configuration. The bin size is AT = 1/10. An excellent agreement is 
found with the limit law (^^, with p = 1/2, hence a = e'^/^/A = 0.151633. 

2.4 Spin correlations 

The time- dependent spin correlation function reads 

Cnit) = {aocrn)t = ((1 - 2ro)(l - 2r„))t = 1 - 4ci(t) + 4rfi,„_i,i(t), (2.28) 

where we have introduced the one-cluster function c„(t) and the two-cluster function 
dm,k,nit), defined as 

Cnit) = (ri . . . Tn)t = n(o_^), (2.29) 

n 

dm,k,n{t) = in. . . TraTm+k+1 ■ ■ ■ Tm+k+n)t = Vtio-'-O ■■ ■ O - ■ - o). (2.30) 
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X 0.2 - 




Figure 2: Distribution of the jamming time T/v of CIC samples of = 1000 spins, against X]^ 
defined in (2.26). Histogram: numerical data {p = 1/2). Thick line: limit Gumbel law ( |2.27| ). 



The one-cluster function is the probabiUty that the sites 1, . . . , n belong to a cluster 

of at least n consecutive down spins. It is therefore directly related to Pn(^) defined 
in (I2I5I) by 

Pn{t) = Cn{t) - 2c„+i(t) + C„+2(t). (2.31) 

Equation ( 2.19 ) then yields 

Cn{t) = A{t) z{ty\ (2.32) 

with (|2:20| ) and 

ait) t ,,-t 



e exp{p{e - 1)). 



(2.33) 



(1 — pe *)^ 

This result can be alternatively recovered by deriving rate equations for the one-cluster 
function itself. As a consequence of ( |2.2| ), each variable r2,...,r„_i entering the def- 
inition ( [^.29| ) can flip from 1 to with rate c„(t) per unit time, while the rate for 
the first variable ti is ac„(t) + (1 — a)c„+i(t), and the rate for the last variable r„ is 
acn+i{t) + (1 — a)c„(t). Gathering all the contributions, we obtain the equation 



dcnjt) 
dt 



■{n - l)c„,(t) - c„+i(t). 



(2.34) 



whose solution coincides with (|2.32|), (BT 
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A similar analysis yields the following rate equations for the two-cluster function: 



dt 



[m + n- 2)dm,k,n{t) 



(2.35) 



0'{dm+l,k,n{t) + dm,k-l,n+l{t)) — (1 — a){dm+l,k~l,n{t) + dm,k,n+l it)) 



for m,k,n > 1, with initial conditions dm,k,n{0) = and boundary values dm,o,n(t) 

Cm+n{t). The rate equations ( p. 351 ) are solved by the Ansatz 



dm,kAt) = Mt) ^t)"^^"", 

provided the amplitudes Bk{t) obey 



dt 



{2- z(t))Bk{t)- z{t)Bk-i{t) 



(2.36) 



(2.37) 



for /c > 1, irrespective of the value of a, with the initial condition Bk{0) = 1, and the 
boundary value Bo{t) = A{t) (|CT| ). 

In order to solve (|2.37|) , we introduce the generating series 



Bix,t) = J2Bk{t) 



X 



(2.38) 



k>l 



which obeys the differential equation 
dB{x,t) 



dt 



(2 - (x + l)z{t))B{x, t) - xz{t)A{t), 



considering x as a parameter, with initial condition i3(x, 0) = x/(l — x). This equation 
can be solved by 'varying the constant': 



B{x,t) =e2*exp(p(e~* - 1)) 



— H ^ — I exp(pa;(e * - 1)) ^ — e * 

px I — X J px 



(2.39) 



Inserting ( p.32| ) and (|2.36|) into ( |2.28| ), we obtain 

Cnit) = 1 - 4pexp(p(e-* - 1)) + Ve-2*5„_i(t). 



Finally, expanding ( |2.39|) and using ( p.22| ), we obtain the expression of the connected 
spin correlation function: 

cr"(<) = c„m - M{ty- 



m>n+l 

As a consequence, in the blocked states, the correlation function reads 
Cr"(oo) = C„(oo)-M(oo)2 = 4pe-M(l-p)t^-p ^ ^"^^ 



ml 



m>n+l 



(2.40) 
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The first term in the above expressions is the leading one, implying that the connected 
correlation function has a factorial asymptotic decay of the form {p{l — e~*))"/n!, mod- 
ulated by an oscillating sign, for any value of p and any time t. This super-exponential 
fall-off is a characteristic feature of irreversible processes such as RSA [^. This be- 
havior is entirely missed by the a priori approach (|2.11|) , (p.l2|) , where correlations fall 
off exponentially, as they generically do in equilibrium systems. Figure ^ shows a loga- 
rithmic plot of (— !)"■ C™"'^(cxo) against n, for a random initial configuration [p = 1/2), 
together with both predictions of the a priori approach, i.e., the full ensemble (p.l2|) , and 
the restricted ensemble ( |2.11| ) where the exact magnetization ( p.24|) is imposed. Both 
predictions appear as straight lines on the plot. The exact value of C5^°°'^(oo) is correctly 
reproduced in the restricted a priori ensemble. 




Figure 3: Connected spin correlation function in the final states of the CIC. Full (open) sym- 
bols show positive (negative) correlations. Circles and full line: logarithm of (—1)" times 
the connected correlation C^°'^'^(oo) ( |2.40 ) for p = 1/2, against n. Squares and dashed line: 



prediction ( 2.11 ), ( 2.24 ) of the restricted a priori ensemble. Triangles and dashed line: predic- 
tion ( p. 12 ) of the full a priori ensemble. 



2.5 Distribution of final magnetization and dynamical entropy 

We finally determine the full distribution of the number of spin flips and of the final 
magnetization, for a given finite sample. This problem was tackled long ago by a some- 
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what similar approach in the case of dimer deposition, without consideration of the 
dynamical entropy, though. 



Single active cluster 

We consider first the case of a single active cluster of size £ > 2, with free boundary 
conditions, with all spins being initially down. In the language of deposition, this corre- 
sponds to an initially empty cluster. We are interested in the distribution of the number 
of spin flips (i.e., deposited particles) during the history of this cluster, until it reaches 
a blocked configuration. The final magnetization Mi of the cluster is such that 

= ^a„(oo) = 2z/,-£. 

n=l 

Because of the irreversible character of the dynamics, every spin flip replaces the 
cluster where it takes place by one or two smaller clusters, according to ( p.l7|) . In the 
generic case, two clusters of lengths di and £2 are generated, and their subsequent histories 
are independent. We have therefore 

z/^ = l + z/£, +z/£,, (2.41) 

where v^^ and are independent random variables, whose distribution is to be deter- 
mined, while = n — 1 and £2 = i — n, with the breaking point n being uniform in the 
range 2 < n < i — 1. For n = 1 ot n = i, only one cluster is generated, and ( |2.41| ) is 
changed accordingly. Finally, we set 1/0 = 0, which necessarily holds, and z/i = 0, which 
contains the gist of the kinetic constraint in the CIC model. Let 

0,(A) = (e^'^O (2.42) 

be the characteristic function of the distribution of Ui. Equation ( p.41|) implies 

e-i 

(^-l)0,(A) = e"^0,(A)0,_fe_i(A) 
fc=i 

for i > 2, with (polX) = 0i(A) = 1. These quadratic recursion relations can be solved by 
introducing the generating series 

$(x,A) = ^0^(A)a;^ (2.43) 



which obeys 
with $(x, A) = l + x + -- -asx— >0 



^d$(x^ = ($(x. A) - 1)(1 + xe^<!>(x, A)), (2.44) 
ax 
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The quadratic differential equation ( |2.44|) has an obvious solution A) = 1. Setting 
A) = 1 + l/u{x, A), we obtain a linear equation 

d-ufx, A) ,^ \n / ,N \ 
X — + (1 + xe^)u{x, A) = -xe^, 
ax 

which can be solved by 'varying the constant'. This yields 

^^ exp(xe^) + e'*^ - 1 .r^ 

* ^) = T. T^^ 7 — \ 7- 2.45 

^ ' (1 - xe^) exp(a;e^) + e^ - 1 ^ ' 

This expression formally contains the distribution of the number z/^ of spin flips. 

First, by expanding ( p.45| ) around A = 0, we obtain generating series for the successive 
moments of v^. The first of these series, 

can be inverted explicitly, yielding 

= (l-e-^K-e-i + (£+l) E ^2.46) 

The mean number of spin flips therefore grows linearly with the cluster size with a 
coefficient Fi = 1 - e"^ = 0.632120, in agreement with ( p.25|) . The constant — e ^ can 
be viewed as the contribution of the free ends of the cluster, while the last term falls 
off factorially, with an oscillating sign, as (— 1)^/(£ + 1)!, just like the connected spin 
correlation ( p.40| ) for p = 1. 

Then, by inverting the generating series ( |2.43D , we obtain the exponential estimate 

0,(A)~e^^W, (2.47) 

where Xc(A) = exp(— F(A)) is the zero of the denominator of (|2.45|) . As a consequence, 
all the cumulants of vn grow linearly with the size I of the cluster, as 

with 

We recover the above result {yi) ~ F\i^ whereas Varz/£ «i ^2^, with F2 = 3e ^ — e ^ = 
0.038126. The bulk of the distribution of vi is therefore a Gaussian of the form 

P(.,) ^ (2vrF2£)-/2 exp [- ^"''Jf \ ■ (2-48) 
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In order to investigate the tails of the distribution of ui for (. large, we set 

^~ I ~ 2 ■ 
An inverse Laplace transform of ( p.42|) , using ( |2.471) , yields 



P(z/,) ~ j 



^ ^t(F(\)-\i) 

27ri 



Evaluating this integral by the saddle-point method, we obtain an exponential estimate 
similar to ( |2.8| ): 

P(z/,)~exp(-£S(0), (2.49) 
where the functions F{X) and S(,^) are related to each other by a Legendre transform: 

dy d P 

F(A) + E«, = A^. A=-. « = -. 

The function S(^) is the large-deviation function (or entropy function) of the quantity z/^. 
This is a positive, convex function of ^, which vanishes quadratically around (^) = Fi, as 



2F2 

in agreement with the Gaussian law ( p.48|) . 

Coming back to the language of the magnetization M, the above functions have the 
following parametric form, in terms of 2 = Xc{X)e^: 

F = ln(l - (1 - z)e') - \nz, X = ln(l - (1 - z)e'), 

E^ln.-'-",r'°'ln(l-(l-.)e'). M = 1 - 2 ' " ".r'^' ■ 

z^e"" z e 

Uncorrelated initial state 

We now investigate the distribution of the final magnetization for a finite chain of 
spins, with an initial state of the form ( ^^.131 ). This magnetization is given by 

N 

NMn = 2u + NMNiO) = 2u+Y,crn{0), (2.51) 

n=l 

where Mjv(O) is the initial magnetization and u is the number of spin flips during the 
history of the system. The final magnetization is therefore random in two respects, 
as it depends both on the initial spin configuration and on the numbers of spin flips 
during the history of each cluster. 

We again introduce the characteristic function 

ijNiX) = {exp{XNMN)) = (exp (^2Az/ + A J] a„(0)^ ^ (2.52) 
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as well as the generating series 



vl>(x,A)= E^iv(A) 



X 



N 



N>1 

The brackets in the right-hand side of ( |2.52| ) involve: 

(i) averaging over stochastic histories with a fixed initial configuration, 

(ii) averaging over the distribution (|2.13| ) of initial configurations. 

The outcome after (i) is that the right-hand side of ( p. 52 ) is a multiplicative cluster 
quantity of the type investigated in Appendix A, where the contributions of clusters of 
up and down spins read 

/i = e^^, ^i = e-^V(2A). (2.53) 
Step (ii) can now be performed. The generating series corresponding to (|2.53|) are 

xe'^ 
1 — xe 

Using ( |2.45| ) and ( [A.l[ ), we obtain 



/(^) = T^^' 9{x) = $(xe-\ 2A) - 1. (2.54) 



xe'^ (exp(pxe^) + (1 — p)(e'^^ — 1] 

^(x, A) = ^ 1 ^ — -I . (2.55) 

(1 — xe^) exp{pxe^) + (1 — (1 — p)xe^){e'^^ — 1) 

This expression provides the distribution of the final magnetization, for any system 
size and any value of the parameter p characterizing the initial state. 

By expanding ( |2.55| ) around A = 0, we obtain generating series for the moments of 
NMj^f. The first of these series yields an expression similar to ( |2.46D for N{M]\f), with a 
leading term, linear in A^, in agreement with the expression (|2.23| ) of M(oo), a constant 
boundary term, and an oscillating, factorially decaying correction. Similarly, we find 

A^VarM^4pe-P((2/_p + 2)e-P- 1). 

For a random initial configuration {p = 1/2), we have therefore A^VarM ^ 4e~^ — 
2e-i/2 = 0.258456. 

Table 1 provides a comparison between exact dynamical results for a random ini- 
tial state {p = 1/2) and prediction of the full a priori ensemble, concerning the main 
characteristics (mean value and scaled variance) of the final energy of the three models 
considered in this work. 

The tails of the distribution of are again described by an exponential estimate of 
the form (1^ ): 

P(M^) ~exp(-A^Sp(M;v)), 
where the large-deviation function Sp(M) reads, in parametric form: 

{l-{l-p)z){l-{l-p)z-{l-z)eP') 1- {l-p)z- {1- z)eP' 

— 111 Z — ~, ; ; — ; ill ■ 



^ pz^(2 — p — (1 — p)z)eP^ 1 — (1 — p)z ' 

pz^{2-p - {l-p)z)eP'^ ■ ^ ■ ^ 
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Model 


E 

dynamical 


E* 
a priori 


AT Var E 

J. y V CAii- -I / 

dynamical 


AT Var E 
a priori 


CIC 
Glauber 
Kawasaki 


-0.393469 
-0.632121 
-0.274087 


-0.447214 
-0.447214 
-0.236840 


0.258456 
0.406006 
0.459839 


0.357771 
0.357771 
0.527638 



Table 1: Mean value and limit scaled variance of the final energy: comparison of the exact 
dynamical results for a random initial condition (p = 1/2) with the prediction of the full a 
priori ensemble. 



This function has finite limits 

E,(0) = -1 In E,(l) = - ln(l - p), 

at the minimum magnetization M = 0, corresponding to z = 0, and at the ground-state 
magnetization M = 1, corresponding to z = 1/(1 —p). Furthermore, the result ( p.50| ) is 
recovered by setting p = 1 in (|2.56|), as it should be. 



Figure ^ shows a plot of the dynamical entropy, defined as being the large- deviation 
function Sp(M) of ( |2.56| ) for a random initial configuration, i.e., p = 1/2, against the 
magnetization M of the final state. The endpoint values read Si/2(0) = ln(8/3)/2 = 
0.490415 and Si/2(1) = ln2 = 0.693147. The prediction ( |2.9|) of the a priori approach 
is plotted for comparison. The functions Si/2(M) and Zlap(M) respectively vanish for 
M{oo) (|2.24| ) and M* (|2.7| ). These numbers are listed in Table 1, together with the 
corresponding limit scaled variances. 



3 Constrained Glauber dynamics 
3.1 Definition of the model 

We now consider a ferromagnetic Ising chain with Glauber dynamics (non-conserved 
order parameter) in the presence of kinetic constraints. The Hamiltonian of the chain, 
with unit exchange constant, reads 

= - ^ (^n<7n+l = - ^Sn, (3.1) 

n n 

where we have introduced the energy (bond) variables s„ = cr„(T„_|_i. 

We consider single spin-flip (Glauber) dynamics, assuming that the flipping rate only 
depends on the energy difference between the configurations after and before the proposed 
move, i.e., 

W{an -an) = Wsm 



19 




Figure 4: Full line: plot of the dynamical entropy of the CIC, given by T,i/2{M) ( |2.56| ), against 
magnetization M. Dashed line: prediction ( |2.5D , ( |2.9D of the a priori approach. 



with 



6n = 2(a„_i + (T„+i)(T„ = 2(s„_i + Sn) e {-4, 0, 4}. 



The requirement that the dynamics obeys detailed balance with respect to the Hamilto- 
nian at temperature T = yields a single condition: 



W4_ 



-4/3 



Choosing time units such that VV_4 = 1, we have W4 = e~^'^. We restrict ourselves to 
zero-temperature dynamics, so that W4 = 0. The rate Wq, corresponding to diffusive 
rearrangements at constant energy, remains a free parameter. The zero-temperature 
limits of the Metropolis and heat-bath rules respectively correspond to Wo = 1 and 
Wo = 1/2. Here we choose 

Wo = 0, (3.2) 
so that only spin flips which lower the energy are allowed. The condition (|3.2| ) defines the 



constrained Glauber dynamics already considered in ||I8| , p^ . The possible spin moves 
are flips of isolated spins: 



+ - + ^ + + +. 



(3.3) 
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Each move suppresses two consecutive unsatisfied bonds: s„_i = s„ = — 1 — s„,_i = 
Sn = +1. The system eventually reaches a blocked state, where there is no isolated spin. 
Equivalently, each unsatisfied bond (or domain wall) is isolated. Our aim is again to 
provide a statistical description of the blocked states reached in this way. 

We recast the problem in terms of deposition, where empty sites represent unsatisfied 
bonds, while occupied sites represent satisfied bonds: 



The moves (]3.cj|) read 



oo 



so that the dynamics is equivalent to the RSA of dimers, considered long ago ^ . 

The blocked states are the spin configurations where unsatisfied bonds are isolated. 
These blocked configurations are therefore formally equivalent to those of the CIC of sec- 
tion 2, up to the replacement of the spins o"„ by the energy variables s„. The Hamiltoni- 
ans and (|3.1|) are also equivalent, up to the replacement cr„ — > s„. As a consequence. 



the entropy S'ap(-E) of the a priori ensemble at fixed energy E is still given by (|2.5|) , up 
to the replacement of M by —E. 



3.2 Dynamics of cluster densities and energy 

We again consider an initial state similar to (|2.13| ), with (To(0) = ±1 at random, while 
each energy variable is drawn from the binary distribution 

f s„(0) = -1 (o) with prob. p, 

\s„(0) = +l (•) with prob. 1 - p. ^'^ 

The parameter p is related to the initial energy E{0) = — 1 + 2p, and (|2.14| ) still holds. 

The dynamics of the cluster densities and energy can be investigated by the method 
of section 2.3. The densities pe(t) of clusters of exactly £ consecutive unsatisfied bonds 
(empty sites) obey linear equations similar to ( |2.18| ): 

^ = -{i-l)pe{t) + 2 E p,it) 

k>e+2 

for i > 1, with Pi{0) = (1 — pYp''', and the energy reads 

E{t) = -l + 2j2ipe{t). (3.6) 

£>1 



The Ansatz (|2.19|) again holds, yielding the solution 

p^{t) = (1 -pe-*)2exp(2p(e-* - 1)) p'e-'^'-'^' 
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and 

E{t) = -1 + 2j9exp(2p(e-* - 1)). 
Again, only inactive clusters of length i = 1 survive in the blocked states: 

Pi(oo) = pe"^^, 

so that 

^(oo) = -1 + 2pe-^P. (3.7) 
This result [ IS , 1^ was shown in Figure |l[ 



For an initial state close to the ferromagnetic ground-state {E{0) —>■ —1, i.e., p —>■ 0), 
the behavior E{oo) ~ E{0) — is easily explained in terms of clusters of two unsatisfied 
bonds. The energy of blocked states then increases monotonically against p, up to the 
maximum value 

^(oo)p=i/2 = -1 + e"^ = -0.632121, 

corresponding to a random initial configuration {p = 1/2, i.e., E{0) = 0), and then 
decreases monotonically against p, down to the value 

E(oo)p=i = -1 + 2e-2 = -0.729329, 

corresponding to the antiferromagnetically ordered initial state {p = 1). 

3.3 Spin and energy correlations 

In the present context, it is natural to consider the spin (site) and energy (bond) corre- 
lation functions 

Cn{t) = ((ToO"n)t, r„(t) = {SoSn)t = {(^0(^lCrn<7n+l) t- 

The energy correlation function r„(t) can be evaluated analytically, using the method 
of section 2.4. We introduce variables r„ = (1 — s„)/2, and consider the one-cluster 
function c„(t) and the two-cluster function dm±,n(t), defined in ( p. 291) and (|2.30|) . These 
functions obey rate equations similar to (|2.34|) and ( p.35|) : 



"^"""^^^ =-(n-l)c„(t)-2c„+i(t). 



dt 

ddrrijk,?! (^) 



{m + n- 2)dm,k,n{t) 



dt 

~ dm+l,k,n(t) ~ dm+l^k-l,n(t) " C^m,fc-l,n+l (^) ~ dm,k,n+lif)- 

After some algebra we are left with the following expression for the connected energy 
correlation function r™'^°(oo): 

r--(oo) = r„(oo)-i?(oo)^ = 2pe-2M(l-2p)tM!-2p , (3.8) 
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which closely resembles ( 2.40 ). For p 7^ 1/2, the first term is leading, hence r™'^°(oo) ~ 
{—2pY/n\. Yoi p = 1/2, r™°'^(oo) ~ (— l)"/(n + 1)!. Figure || shows a logarithmic plot of 
both correlation functions, against n, for a random initial configuration {p = 1/2). The 
circles show (— l)'^ r™°'^((X)), as given by the analytical result ( |3.8| ). The triangles show 
the full spin correlation function C„(cxo), measured in a numerical simulation. For each 
sample, starting in a random initial configuration, the constrained dynamics is run until 
a blocked state is reached. The correlation function C„(cxd) is found to be positive and to 
decay monotonically to zero as a function of the distance n. The data shown correspond 
to a total of lO^'' blocked spins. Both correlations are observed to fall off as l/(ri + 1)!. 




Figure 5: Spin and energy correlation function in the blocked states of the ferromagnetic chain 
with constrained Glauber dynamics. Full (open) symbols show positive (negative) correlations. 
Circles and full line: logarithm of (—1)" times the connected energy correlation r5f°'^(oo) ( |3.8D 
for p = 1/2, against n. Triangles and full line: logarithm of the full spin correlation function 
Cn(oo), against n, measured in a numerical simulation. Dashed line: logarithm of asymptotic 
behavior l/(n + 1)!, up to a multiplicative constant, meant as a guide to the eye. 



3.4 Distribution of final energy and dynamical entropy 

We now investigate the distribution of the number of spin flips and of the final energy, 
using the method of section 2.5. 

We consider first the case of a single cluster of size i > 2, whose initial configuration 
is antiferromagnetically ordered, i.e., made of unsatisfied bonds. Let be the number 
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of spin flips during the history of this cluster. Equation (|2.41|) is replaced by 

ft = 2 + Pii^ + vi^ , 

where i\ = n — \ and i2 = i — n — \^ and the breaking point n is uniform in the range 
\ < n < I — \. Hence equation ( p.44|) for the generating series A) is replaced by 



d'l'(x,A) 
dx 



(xe^$(x,A))%<l>(a;,A) -1. (3.9) 



This is a Riccati equation, which can be solved by linearization [^]. Setting 

1 ^ X du{x, A) 



A) u{x, A) dx 

yields 



d u{x, A) 



e^^ui^x, A), 



dx"^ 

so that a basis of solutions reads exp(±xe'*'). We thus obtain the closed- form expression 

(e^ + 1) exp(2xe^) + e^ - 1 



$(x,A) 



(e^ + 1) (1 - xe^) exp{2xe^) + (e^ - 1) (1 + xe^) ' 



We now consider the final total energy NEn of a finite chain of spins, with the 
initial state (|3.5|) . Equation (|2.51|) is replaced by 

NEn = NEn{0)-2u, 

where En{0) is the initial energy and u is the number of spin flips. Equations ( ^.541 ) 
and (|2.55| ) are replaced by 

fix) = \ , g{x) = ^xe\ -2A) - 1 
1 — xe 

and 

X ((e^^ + 1) exp(2pxe-^) + {2p - l){e'^^ - 1)) 
~ (e^^ + l)(e^ - x) exp{2pxe-^) + (-e^ + (1 - 2p)x){e^^ - 1)' 
This last expression contains the distribution of the final energy of a system of size N, 
as a function of the parameter p characterizing the initial state. In particular, the mean 
energy is found to agree with the expression ( p.7[ ) of E{oo), while its scaled variance 
reads 

NVarE ^ Ap^ip'^ -p+ l)e-^P. 

For a random initial configuration {p = 1/2), we have therefore NVaiE ^ 3e~^ = 
0.406006. 

The tails of the distribution of En are again given by an estimate similar to ( |2.8| ): 

P(E^) ~exp(-iVSp(E^)), 
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where the large-deviation function reads, in parametric form: 



(1 + {2p - l)zf -{z- l)^e^P' l + {2p- l)z + (1 - z)e'^P' 



-1 + 



2p^2g2p^(2(l - p) + (2p - l)z) 
This function has finite hmits 



(3.10) 



Sp(-l) =lnz,(p), Ep(0) = --ln(p(l-p)), 



at the ground-state energy E = —1, corresponding to 2; = Zc{p), with 

1 + (2p - l)z, + (1 - z,)e^P'^ = 0, 



(3.11) 



and at the maximum energy E = 0, corresponding to 2; — > 0. Figure ^ shows a plot of the 
dynamical entropy Si/2(-E), as given by ( p.lO| ), against the energy E of the final state. 
The prediction of the a priori approach is plotted for comparison. The endpoint values 
of the dynamical entropy are Si/2(— 1) = ln2;c(l/2) = 0.245660 and Si/2(0) = ln2 = 
0.693147. 




Figure 6: Full line: plot of the dynamical entropy of the ferromagnetic chain with con- 
strained Glauber dynamics, given by Y^i/2iE) ( |3.10 ), against energy E. Dashed line: pre- 
diction (^), ( ^ ) of the a priori approach. 
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4 Constrained Kawasaki dynamics 



4.1 Definition of the model 

We finally investigate a ferromagnetic Ising chain with conserved dynamics at zero tem- 
perature, in the presence of kinetic constraints. 

Consider the ferromagnetic chain with Hamiltonian ( |3.1|) , but now with Kawasaki 
dynamics, where only pairs of opposite spins (s^ = anCn+i = —1) may be flipped, so 
that the magnetization is locally conserved. The flipping rates are again assumed to 
depend only on the energy difference involved, which now reads 

6H = 2(a-„_ia-„ + crn+icrn+2) = 2(s„_i + Sn+i) G {-4, 0, 4}. 
We make the choice (|3.2|), defining thus a model with constrained Kawasaki dynamics, 



already considered in |2y, |2T|, ^ . The possible spin moves are 

- + - + ^--++, +- + > + + --. (4.1) 

Each move suppresses two unsatisfied bonds: Sn-i = Sn+i = —1 ^ = = +1. 

The system eventually reaches a blocked state, where the spin patterns H 1 — and 

— I h are absent. Equivalently, there are at most two consecutive unsatisfied bonds. 

Considering unsatisfied bonds as empty and satisfied bonds as occupied, as in ( |3.4|) , 
the moves ( [4.1D read 

000 — !► • O •, 

SO that the dynamics is equivalent to the random deposition of hollow trimers, a case of 
RSA that seems not to have been studied so far. 



4.2 A priori statistics 

We consider the restricted a priori ensemble of blocked configurations with energy E. 
The entropy S'ap(-E') of this ensemble can again be evaluated by means of the transfer- 
matrix formalism. The partition functions of a finite chain of spins, now labeled by 
the prescribed values of its last two bonds (o or •), obey the recursion 

where the 4x4 transfer matrix 



T = 



( e^ 







o\ 


e-^ 





e-^ 













e^ 


I 
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has a reducible characteristic polynomial, AP3(A), with -P3(A) = — e'^A^ — A — e~^. 
Thermodynamic quantities are still given by ( |2.10| ) in terms of the largest root A^ of P3, 
i.e., in parametric form: 



hence 

l-z^ _ {z^ + z + 1) ln(2:^ + 2: + 1) - z{2z + 1) In z 

The entropy of the full a priori ensemble of blocked states, irrespective of their energy, 
is equal to the maximum value of the entropy S'ap(-E), corresponding to /? = 0, where we 
have — — z — 1 = hence z = Zq = 1.839287 and 

S* = In^o = 0.609378, E* = ^ = -0.236840. 

''P z^ + 2zo + 3 



The difference Sap(-E') = S'*p — S'ap(-E'), introduced in (p.9[), to be plotted in Figure 
vanishes quadratically as 



^2 „ _ (zg + 2zo + 3)^ 
84 (zo' + 4^0 + 1) 



so that A^VarE ^ l/(2c) = 0.527638. 

4.3 Dynamics of cluster densities and energy 

We again consider an initial state of the form ( p.5| ). The densities pe{t) obey linear 
equations similar to ( P^.18| ): 

%^ = -(£-2)Mt) + 2 E Mt) 
^'^ k>e+3 



for i >2, with ^^(O) = (1 — p)'^p'''. The Ansatz (|2.19| ) yields the solution 



p,(t) = (1 - pe-'Y exp (2p(e-* - 1) + p\e'^' - 1)) p^e'^'-^^'. (4.3) 
The dynamical equation for £ = 1, 

^=p,{t) + Y.kpdt), 

^'^ fc>4 

is special, because any move generates an isolated unsatisfied bond. Using (^73|) , we get 

pi{t) =p + p\pe-'-2)exp(2p{e-'-l)+p\e-^'-l))-2p^e-'^^+P^' J^^" dy. (4.4) 
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Equation (|3.6|) then yields 

r-l+p 
'1+pc- 

Only inactive clusters of one or two unsatisfied bonds survive in the final states: 



E{t) = -1 + 2p - Ve-(^+P)' f^^" e^' dy. 



pi(oo) =p- 2p2e-W _ 2p2e-(i+p)' J^^^"' dy, ^2(00) = /e-'^"*'', 

so that 

^(00) = -1 + 2p - Ve-(i+P)' ^^^"^ e^' dy. (4.5) 

This result was shown in Figure |I]. 

For an initial state close to the ferromagnetic ground-state (-E'(O) — > —1, i.e., p 0), 
the behavior £'(00) ^ E{0)—4p^ is easily explained in terms of clusters of three unsatisfied 
bonds. The energy of blocked states then increases monotonically against p, to the 
maximum value ^ 

E{oo)p=i = 1 - 4e-^^ e^' dy = -0.098204, 

corresponding to the antiferromagnetically ordered initial state {p = 1). For a random 
initial configuration {p = 1/2), we have 

^(oo)p=i/2 = -e-^/^y e^ dy = -0.274087. 



The last two results are already in EQ, 22 



4.4 Spin and energy correlations 

The energy correlation function r„(oo) can be evaluated analytically, similarly to sec- 
tions 2.4 and 3.3. In the present situation final results are, however, less explicit. 

We consider the one-cluster function Cnit) and the two-cluster function dm,k,n(t), de- 
fined in ( |2.29|) , ( p.30|) . The rate equations obeyed by these functions, and the way to 



solve them, are similar to ( p.34| ), ( p.35| ) 
The one-cluster function obeys 



= -{n - 2)c„(t) - 2c„+i(t) - 2c„+2(t) {n > 2) 

dci(t) 



-2C3(t), 



hence 



dt 



Cn(t) = exp {2p{e-' - 1) + p\e''' - l)j p^e"^""'^* (n > 2) 
ci(t) =p- 2p2e-(i+P)' Z'^" e^' dy, 

Jl+pe-' 
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so that 

ci(oo) =p- 2p2e-(i+P)' dy, 02(00) = p^-^^-P', (4.6) 

in agreement with ( P?^ ), P?^), and P^). 
The two-cluster function obeys 

~ d m+l ,k,n{i) ~ drn+l,k-l,n{t) " dm,k-l, n+l (t) ~ dm,k, n+l {t) (4-7) 
~ d m+2 ,k,n(t) ~ d m+2 ,k-2,nit) — dm,k-2, n+2 {t) — dm,k, n+2 (t) . 

Besides the conventions of section 2.4, only some of the terms are present in the right- 
hand side for either m or n = 1 or 2, namely those for which the underlined index is 
greater than 2. Furthermore, for A; = 1 the sum dm+2,k-2,n(t) + dm,k-2,n+2it) is replaced 



by Cm+n+iit). As in previous cases, we look for a solution to (^l7| ) of the form 

dm,kAt) = Ak{t)z{tr+'' {m,n>2), 
dmAiit) = Mt)4tr {rn>2), 
di,k,i{t) = Dk{t), 

where z{t) has been introduced in (|2.20|) . The procedure then consists in introducing 
generating series ^(x,t), i3(x, t), V{x,t), similar to (|2.38|) , writing differential equations 
obeyed by these functions, and solving the latter equations. This requires some lengthy 
and tedious algebra. 

The energy correlation function in the blocked states is still given by (|2.28|) in 
terms of Ci(oo) and (ii^„_i^i(oo), so that the function of most interest is V{x,oo), for 
which we are left with the expression 

V{x, oo) = — p\x + 2)e-(^+P)' [^'^^ e^' dy 

1 — X Jl 



V / o, 2 , 1^ / (x + I + p{x^ + 1)Y 



2{x' + 1) exp ^ ' " / e^ di/ (4. 



1 - X ^ ' ' \ 2{x^ + 1) 



/ dy 

J Jm 



+ — V 2 a; + 1 exp oi 2 ^ i\ / ^0^' « / _ dy, 

X^ ^ \ 2[X'^ +1) J Jo Jb(0) 



with the notations 



R{x, u) = (1 — X + px{l — X) + 2p X ) 

1 — X 

2 



X exp \^p{x + l)(n - 1) + y (x^ + 1)(m^ - 1) 
(1 — + px{l — 2x — x^)u) 

X exp \ p{l - x){u - 1) + y (1 - x'^){u^ - 1) 
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+ {x" + 1) / dy 

J l+pxu 

X exp {l+p{x-l + {x + l)u) + y (a;^ - 1 + (x^ + 



and 



, X + 1 + l)u 

b[u) = — 



2(x^ + l] 



The function V{x, oo) has a simple pole at x = 1, with residue — ci(cxd)^, where ci(oo) 
has been evaluated in ([4.6D, so that the di^k,i{oo) converge to Ci(cxd)^ as the distance k 



becomes infinitely large, as it should. The fall-off of the difference (ii,A:,i(oo) — Ci(oo)^, 
and that of the connected correlation function r™'^°(oo), are related to the behavior of 
V^x, oo) as is large. The result (^78|) leads to the estimate 

V^x, oo) ~ exp{2px — p^x^)^ 

with exponential accuracy. The results summarized in Appendix B then imply 



The result ([4.8|) does not however lead to any useful expression for r„, even for n = 1. 
Figure ^ shows a logarithmic plot of the spin and energy correlation functions against n, 
as measured in a numerical simulation for a random initial condition (p = 1/2). The data 
shown correspond to a total of 5 ■ 10^° spins. Both the full spin correlation C„(oo) and 
the connected energy correlation r™'^°(oo) = rn(oo) — E{ooY are found to agree with 
the asymptotic result (f4.9|) . The absolute value of the data follows the predicted fall-off, 

shown as a dashed line, while the signs roughly follow the predicted pattern (+ H ), 

up to more and more seldom mistakes. 



4.5 Distribution of final energy and dynamical entropy 

We end up with the distribution of the final energy of a finite sample. This analysis will 
follow the lines of section 2.5 and 3.4, the main difference being that ([4.10| ) will have to 
be solved numerically. 

We consider first the case of a single cluster of size i >2, whose initial configuration 
is only made of unsatisfied bonds. Let ui be the number of spin fiips of this cluster. 
Equation ( |2.41D is replaced by 

Z/^ = 2 + + , 

where ii = n — 1 and £2 = i — n — 2, and the breaking point n is uniform in the range 
1 < n < i — 2. Hence equation ( p.44| ) for the generating series A) is replaced by 

^d$(x^ = a;V^$(a;, A)^ + 2$(a;, A) - a; - 2. (4.10) 
dx 
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Figure 7: Spin (site) and energy (bond) correlation functions in the blocked states of the 
ferromagnetic chain with constrained Kawasaki dynamics, measured in a numerical simulation. 
Pull (open) symbols show positive (negative) correlations. Circles and full line: logarithm of the 
absolute value of the connected energy correlation r^°°'^(cxD). Triangles and full line: logarithm 
of the absolute value of the full spin correlation C„(oo). Dashed line: logarithm of asymptotic 
behavior l/(2"(n/2)!) up to a multiplicative constant, meant as a guide to the eye. 



In contrast with ( p.44| ) and ( |3.9|) , we have not been able to solve the Riccati equa- 
tion ( [4.10|) analytically. 

We now consider the final total energy NE^ of a finite chain of N spins, with the 
initial state (|3.5|) . The tails of the distribution of E^ are again given by an estimate 
similar to (|2.8|), 

P(E^) ~exp(-iVSp(E^)), 
where the large-deviation function ^p{E) reads, in parametric form, 

^ = -^dA' ^^ = ^""^-^dA' ^^-''^ 
and Xc{X) = exp(— F(A)) is the real positive solution of 

$(px,e\-2A) = — (4.12) 

(1 -p)Xc 

The function Sp(i?) can be evaluated anal5d;ically in some regimes. Skipping any 
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detail, we mention that the hmits 



Sp(-l) = - ln(l - p), Ep(l/3) = -- ln(p2(l _ p)), 

at the ground-state energy E = —1 and at the maximum energy E = 1/3 can be 
determined exactly. Moreover, the solution to ([4.1CI|) can be expanded as a power series 
^(x. A) = 1/(1 — x) + A<l>i(x) + A^<l'2(a;) + - ■■ We thus recover the mean energy £'(00) ( |4.5|) , 
and obtain the following expression for the scaled energy variance: 

NVaiE ^ 2p(l -j9)(l - V) +4(1 _p)25(p) 

- 4(1 - pY{3 - 4p + Ap^)A{p) - 2(1 - pf{5 - 7p + Ap^)A{pY, 



2p2 2p-l-E{oo) 



with 

{1-py Ji 2(l-j9)2 ' 

B{p) = j^^^-^'^"^'' 1'^' ^'\iy - 2)Aiy - 1) - 2)2 dy. 

Figure || shows a plot of the dynamical entropy Si/2(-E) obtained by solving numer- 
ically ( |4.1CI| ), ( [4.11| ), ( |4.12| ). The prediction ( ^^ ) of the a priori approach is shown 
for comparison. For p = 1/2, we have Si/2(— 1) = Si/2(l/3) = ln2 = 0.693147 and 
A^VarE ^ 0.459839. 



5 Discussion 

We have presented a parallel study of the zero-temperature dynamics of three one- 
dimensional Ising models with kinetic constraints, to which a number of previous studies 
have already been devoted: paramagnetic CIC models [|1^ 0, |1^, [l^ (section 2), fer- 



romagnetic chain with constrained Glauber dynamics [|T8|, |I9| (section 3), ferromagnetic 
chain with constrained Kawasaki dynamics [^, |21|, |22| (section 4). 

The common characteristic feature of these models is that their zero-temperature dy- 
namics is fully irreversible: each spin can flip at most once during its whole history. As 
shown in the present work, these stochastic dynamical systems can be mapped onto mod- 



els of irreversible deposition [^]: CSA of monomers (•) for paramagnetic CIC models, 
RSA of dimers (••) for the ferromagnetic chain with constrained Glauber dynamics, RSA 
of hollow trimers (• o •) for the ferromagnetic chain with constrained Kawasaki dynam- 
ics. This exact mapping onto RSA or CSA models allows the analytical determination of 
many physical quantities. Assuming an uncorrelated initial state with prescribed energy 
(or magnetization), we have obtained exact results for the three models, and compared 
them to the predictions of the a priori approach, testing thus the so-called Edwards 
hypothesis in this particular zero-temperature framework. 
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Figure 8: Full line: plot of the dynamical entropy of the ferromagnetic chain with constrained 
Kawasaki dynamics, against energy E. Data are obtained by solving numerically ( 4.10| ), ( |4.11j ), 
and (4.12), for p = 1/2. Dashed line: prediction (4.2) of the a priori approach. 



We have first shown that the jamming time grows logarithmically with the system size, 
up to finite fluctuations given by extreme- value statistics. The result ( p.26| ), ( p.27|) , estab- 
lished for the CIC, also holds quantitatively for the ferromagnetic chain with constrained 
Glauber and Kawasaki dynamics, with respectively a = p^e"^^ and a = p^e~^P~^^ . 

There is a complete lack of ergodicity in these irreversible models. The mean final 
energy indeed bears a non-trivial dependence on the initial condition, as depicted in 
Figure |l|. For a random initial configuration, the comparison of the exact dynamical 
results for the average and variance of the final energy with the prediction of the a 
priori (Edwards) approach reveals systematic differences, which have either sign, and an 
absolute value ranging up to some 20 percent (see Table 1). 

The two-point spin (site) and energy (bond) correlation function in the final states 
has also been evaluated, either by analytical means or by accurate numerical simulations. 
Connected correlations fall off factorially (see Figures ^ 0); often with an oscillating 
sign. A super-exponential fall-off of correlations is indeed known to be generically obeyed 
in RSA models. Such a feature cannot be reproduced by an a priori ensemble, where 
correlations generically decay exponentially, with a finite correlation length, related to 
the first two eigenvalues of the transfer matrix. 

We have also determined the distribution of the energy of the final states beyond the 
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Gaussian approximation. Such a problem seems to have been tackled only once in the 
RSA literature We thus obtain large-deviation estimates for the exponentially small 



tails of the distribution. The corresponding dynamical entropy depends on the initial en- 
ergy (or magnetization). The comparison of the result for a random initial configuration 
{p = 1/2) with the a priori approach again shows differences at a quantitative level (see 
Figures |, ||, |) . 

The results obtained so far invalidate the Edwards hypothesis in the present situation 
of fully irreversible zero-temperature dynamics. There are indeed systematic differences 
between the exact dynamical expressions and the predictions of the a priori approach, 
and even qualitative discrepancies, such as the super-exponential fall-off of correlations. 

The present work also questions the existence of any simple relationship between the 
landscape of metastable states and the slow dynamics just above the dynamical phase 
transition. Indeed, on the one hand, all the results on the zero-temperature dynamics 
of the CIC are independent of the parameter a, which interpolates between the ACIC 
for a = or a = 1 and the SCIC for a = 1/2. On the other hand, these limiting 
cases are known to have different kinds of slow dynamics in the presence of activated 
processes, at low temperature. For instance, the relaxation time to equilibrium diverges 
exp(2/5) for the SCIC pf, and as Teq ~ exp(/?V(ln2)) for the ACIC 0. 



In spite of its specificity, the present approach may also shed some new light on other 
quantities and/or other situations of interest. One example is the size distribution of 
ordered clusters, which has been recently shown to be a useful tool to test the Edwards 
hypothesis in spin models under tapping |^. In the present context the exact determi- 
nation of the density fe{oo) of clusters of i occupied sites in the final states would require 
a lengthy calculation. However its exponential fall off for a large cluster size is related to 
the ground-state dynamical entropy as fe{oo) ~ exp{—£T,p{E = —1)). We thus obtain 
the simple estimate fiioo) ~ (1 —pY, both for CIC and constrained Kawasaki dynamics, 
expressing that the long ordered clusters in the final state have to be already present in 
the initial state, while the result fe{oo) ~ Zc{p)~^ (see ( p.ll|) ) for constrained Glauber 
dynamics is non-trivial. 

Finally, the present situation of a quench from a disordered initial configuration (in- 
finite initial temperature) can be viewed as the relaxation part of a cycle of random 
tapping with infinitely high intensity. It would also be desirable to extend at least some 
of our results to the more realistic situation of a finite tapping intensity. 
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A Averaging a multiplicative cluster function 

Consider a finite chain of spins cr„ = ±1, numbered n = 1,...,N. Tlie cliain is 
naturally partitioned into clusters of parallel spins. Let M be the number of clusters, 
and Li, L2, . . . , Lm be the lengths of the clusters, with Li + ■ ■ ■ + Lm = N. Assume 
0"! = +1. We have ct„ = +1 for n = 1, ... , Li, then (j„ = — 1 for n = Li + 1, . . . , Li + L2, 
and so on. 

A multiplicative cluster function is a quantity of the form 



Qn 



/Li 9L2 fL;--- if CTl = +1, 
gLifL29L3--- if CTi = -1, 



where each cluster of L up spins brings a factor fi, and each cluster of L down spins 
brings a factor g^. 

Averaging a quantity such as ^at over a state of the form ( p.l3|) amounts to summing 
the contributions of all the partitions of A^ into cluster lengths {L^, k = 1, . . . , M}, with 
the a priori weight 

U,(f r ^^ _ / (1 - P)^^ (1 - P)^' ■ ■ ■ if ^1 = +1' 

^ - \ pL, (1 _ p)L. . . . if ^ 

In order to perform this summation, we introduce the generating series 



N 



L>1 L>1 N>1 

If (Ti = +1, we have 

Qi^) = E fLAa-p)xf' giApxf' fL,{{i~p)xf' ■ ■ ■ 

Lt_,L2,L3,... 

= /((I - P)x) + /((I - p)x)g{px) + /((I - p)x)g{px)f{{l - p)x) H , 

where the successive terms are the contributions of the partitions of the chain into M = 
1, 2, 3, . . . clusters. Adding the contribution of the sector ai = —1, and summing up the 
geometrical series, we obtain the result 

p. X _ /((I - P)x) + gjpx) + 2/((l - p)x)g{px) 

^^""^ I~f{{l-P)x)g{px) ' ^""-'^ 

which interpolates between Q{x) = f{x) at p = and Q{x) = g{x) ai p = 1. 



35 



B Expanding exp{ax — hx^) as a power series 



This appendix is devoted to the power-series expansion 

exp(ax — hx^) = ^ fn{ci, b)x^. 

n>0 

An identification with the generating series of Hermite polynomials p5 



Hn{z)— = exp(22:a; - x^) 

n>0 

leads to 

We are mostly interested in the asymptotic behavior of the coefficients /„(a, b) as n gets 
large, for fixed a and b. The asymptotic expansion of Hermite polynomials [p5| yields 



The above estimate becomes exact for any finite n in the simple situation where a = 0, 
where one has straightforwardly 

/2fc(0,&) = -^, /2fe+i(0,6) =0. 

For generic values of a, the signs of the coefficients fn{(i,b) are given by the cosine 

function in ( |B.1| ). They oscillate according to the four-periodic pattern (+H ), except 

for 'mistakes' which take place more and more seldomly, for n ^ fc^/x, with 



fb 

2^ 



For a > 0, mistakes are isolated + or — signs. For a < 0, they consist of three consecutive 
+ + + or signs. 
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